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On a Method of Finding the Conductivity for Heat 

By Professor C. ISTiven, F.K.S. 

(Eeceived December 5, 1904, — Eead January 26, 1905.) 

Introduction. 

This paper contains an account of a method for finding the conductivity for heat, 
especially in bad conductors. The substance is supposed given in the form of a 
cylinder, or rather of two half-cylinders pressed together and heated by a current 
passing through a wire along the axis. When the steady state has been reached, 
the conductivity is given in terms of the difference of temperatures at known 
distances from the axis, and the heat supplied by the wire. The latter can be 
found when we know the current passing through it and the difference of potential 
at its ends. 

The first part of the paper contains a detailed account of the methods employed 
for finding the difference of temperature, and a description of the apparatus used. 
The results of some experiments made with it are also given, and compared 
with those found by other observers. 

The second part of the paper contains a solution of the mathematical problem 
of the diffusion of heat in an infinite solid from a line at which it is supplied at a 
constant rate, and the solution of some other allied questions. 

One result of the investigation suggests a method of finding the diffusivity 
directly, when the substance is of sufficiently great extent. 

1. If heat be supplied uniformly and continuously to a solid within a 
given space at the rate of H units per second, and if a surface S be drawn 
round this space, then when the steady state has been attained and the 
temperature at each point is constant in time, 

J an 

the integral being taken over the closed surface, h being the conductivity, 
and cWjdn the normal variation of 6 outwards. In particular, if the solid 
be in the form of a circular cylinder, and the heat be supplied by an electric 
current passing through a uniform straight wire coinciding with its axis, 
the equation takes the form 

IL = —k— . 277-r, 

dr 

the enclosing surface being a coaxial cylinder, and B^ the heat supplied per 
centimetre per second. 



1904.] Finding the Conductivity for Heat. 35 

If we integrate this equation, we find for the difference of temperatures 
#!— #2 at distances a, b, from the axis, the expression 

H x log, h/a = 2tt (0 l - 8 ) k 

Upon this equation we may base a method of finding h, since the heat 
evolved in the wire per unit length is readily found when we know the 
current i (in amperes) passing through it and its resistance ; or what comes 
to the same thing, the current and the difference of potentials V 1 — V 3 
(in volts) between two points at a known distance apart. These data can 
be got by using, say, a Weston millivoltmeter provided with suitable shunts 
and resistances, and may be obtained with considerable accuracy if the 
instrument be properly calibrated. 

Reducing electrical to heat units the value of the conductivity is given by 

, i(V l —Y 2 )logeb/a 

I being the length of the heating wire, at the ends of which the potentials 
are taken. 

The heating wires used were of platinoid. At first wires of No. 40 
S.W.G-. were used carrying currents of about one ampere, but these were 
afterwards replaced by thicker wires of No. 34, as they were found to give 
too great a heating effect, and caused the wood in which they were placed to 
shrink at the centre, and besides, in addition, were liable to snap when drawn 
tight. The thicker wires gave more satisfactory results. 

2. The temperatures V # 3 , or rather 6 X — 3 , were measured either by the 
change of resistance of platinum wires embedded in the body parallel to 
the axis, or by thermo-electric junctions formed by soldering together thin 
iron and German silver wires placed in the same positions as the platinum 
wires. These fine wires were stretched round two brass pins so as to be 
at distances about 1 cm. and 3 cm. from the heating wire. 

"When solids were tested they were made in the form of two half 
cylinders, one of which was placed in a semi-cylindrical wooden shell, which 
was moved up till the two wires rested on its fiat face. If platinum wires 
were used their resistances were now compared. The current is then passed 
through the heating wire, which is drawn tight round a brass pin and fixed 
to lie parallel to the other two wires. 

The second half cylinder and the other semi-cylindrical wooden shell 
were then placed upon it and squeezed firmly down. If this operation be 
carefully performed it is possible to fix the wires in their positions, though 
the thin platinum wires gave some trouble owing to their great fragility. 

D 2 
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When the apparatus is used for testing the conductivity of powders, the 
lower shell is first filled and the wires placed in position on the upper 
surface of the substance, the current turned on and the heating wire drawn 
tight. The upper shell is now laid above the first, and filled through a 
slot made along the top, the powder being distributed to fill it completely. 

The apparatus is shown generally in fig. 1. Fig. 2 is a horizontal section 

Fig. 1. 




Fig. 2. 




in the plane of junction of the two half cylinders ; AA, BB are the two 
thermometric wires, and CD the heating wire fixed to brass plates. The 
length of the cylinders used was 15 cm., and the diameter 9 cm. The 
heating wire had a length of 25 cm. in the later experiments. 

3. Thermometric Arrangements. — (1) When the difference of temperatures 
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was measured by the change of resistance of platinum, a piece of platinum 

wire was coated with copper and the whole drawn out till the platinum had 

a thickness of about 0*035 mm., that of the copper being about 0*27 mm. 

Any change of resistance was thus appreciably due to the platinum. Two 

lengths from the same piece were taken, laid together, and being connected 

with soft wax, were bared as nearly as possible to the same extent of about 

3 to 4 cm., so that their resistances were very nearly equal. 

If E be the resistance of the inner one, E' that of the outer before 

heating, let 

E : E' = a : yS, 

when compared by the Wheatstone's Bridge method, taking, say, a = 1000 
and /3 a number not far from 1000. After heating, and when the 
temperature has become steady, 

E becomes E x = E [1-f c(^ 1 — O )], 
R'o » R'i == R'o [1 + c (0 S •— O )], 

where c is the temperature coefficient of the resistance of platinum, and Q 
is the temperature of the room. 

Suppose now that E x : E\ = a : ft^ Then 

If the total changes of resistance are only comparatively small fractions of 
the original resistance, we have 

c (0i-0 3 ) = (/3-&)//3i- 

Of course, if c (# 3 —#o) * s a considerable fraction, the resistances must be 
found separately. In the experiments actually made, 0^—Oq amounted to 
about 30° at most, and c = 1/400, so that the assumption here made is 
sufficiently accurate. 

(2) So much difficulty was experienced in the use of these fine platinum 
wires, especially in testing loose powders such as sand or mould, that they 
were ultimately given up, and the difference of temperatures found by the 
thermo-electric junctions already described. The arrangement used was a 
wire of German silver of 0*2 mm. diameter between two iron wires of the 
same diameter, stretched parallel to the axis of the cylinder so that the two 
wires, and, therefore, the junctions, were at the distances apart previously 
chosen, that is to say, 1 cm. and 3 cm. from the heating wire. As the 
resistance of the galvanometer employed was about 800 ohms, the deflection 
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was sensibly proportional to the difference of temperature. By independent 
experiments one scale division represented 0*70° C. 

4. With regard to the choice of the distances a, b from the central wire, it 
is clear that these should not be taken too nearly equal ; for, unless & be a 
moderate multiple of a, any error in the measurement of one of them, say a, 
will introduce a large error in log* b/a. 

Let a f for instance = a{l — l/n), 

log* b/a / = log b/a+l/n, nearly. 

Thus, if a = 1 cm., b = 1*5 cm., and - a == 0*5 mm., 

n 

log b/a = 0-405, 1/^ = 0-05, 

so that the error is about 12 per cent, of the whole. 

By taking b/a = 3 the error lies between 4 and 5 per cent. 

As between the two pairs of distances 1, 6 and 2, 6 cm., we should rather 
choose the latter, as log* 6 is not quite 2 log e 3, and an error of 0*5 mm. in 
2 cm. is only half as important as the same error in 1 cm. 

In the present form of the apparatus I have therefore taken a = 1 cm., 
b = 3 cm. ; and in a larger model, a = 2 cm., &" = .6 cm. 

5. It is desirable to know how long an experiment need be continued to 
get a fair approximation to the steady state. For this purpose it is necessary 
to know how the temperature rises at any point. 

In the second part of this paper the solution is given of the problem of the 
diffusion of heat in an infinite cylinder from a steady source at the axis. 
Taking the result as a guide for the present case, we have 



0i-e. = \ 



fb/*/nt 
a/ */nt 



where # 2 = r^/nt, and r is the distance from the axis and n the diffusivity. 

Assuming, provisionally, that r*/4m,t is large enough to make r^/^nt small, 
we may expand e" r2/hlt } and retain only the first two terms. Thus 

b W—a?' 



#1-02 = «( l0 g 



a Snt 

The first term of the expression gives the temperature of the steady state, 
and if the temperature be 5 per cent, below its final value, 

If a = 1 cm., 6=3 cm., 

t = 18'2/n. 
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In the case of sand, for which n is given, 0*0032 (Everett's e Units and 
Physical Constants ') 

t = 5702, or about 1^ hours. 

To get within 1 per cent., we should have to continue the heating for about 
7 \ hours. 

It may, perhaps, be worth remarking that the diffusivity may be found 
directly from the expression for ddjdt given in the second part of the paper 
and quoted above ; for 



which is a maximum when 



dt t 



t = r 2 /4:n. 



Taking r = 3 cm., the time required is only about one-eighth of that 
previously necessary. 

6. As an illustration of the results which the method gives, I add a list of 
some substances whose conductivities have been found by it. In the case of 
the woods, the flow of heat is across the fibres. In the second column the 
letter E indicates that the difference of temperatures was found by the change 
of resistance of the platinum wires, J that it was found by thermo-electric 
junctions. The distances of the thermometric wires from the central wires 
are given in the third and fourth columns ; in the earlier experiments, b/a 
was taken nearly = 1-5, while in the later ones b/a is about 3. The 
numbers agree fairly well, though, for the reasons given above, I think the 
latter ratio preferable. The fifth column shows the size of the platinoid wire 
used as the source of heat, and the sixth the current passing through it. 

Some determinations by other observers are added ; these are taken from 
Everett's ' Units and Physical Constants/ 

I wish to express my indebtedness to Mr. William Mitchell, M.A., by whom 
many of these determinations were made, for this and other assistance he has 
given me in the course of the work. 
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Paet IL 

In connection with the foregoing method, based on the steady state, it seemed 
desirable to know the mode in which this state is attained at any time after the 
heating is begun. As a preliminary, the problem of the heating of a rod by 
heat supplied uniformly in the middle, was first attacked in Fourier's way by 
treating it as the limit of a ring of infinite radius. The results, though 
complicated in form, may be worth recording, as the conditions approximate 
to those actually met with when one end of a metal rod is put into a hot 
flame. 

The case of a small spherical cavity is added, and may be taken as an 
illustration of the way in which the heat given out by a small morsel of radium 
is propagated through an infinite solid. 

In all the questions treated in this part the body is supposed infinite in the 
directions in which the flow of heat is considered, the ring and thin spherical 
shell excepted. 

Case 1. Thin Circular Ring with Surface Impermeable to Heat — Consider 
first the case of a uniform thin circular ring heated at one spot. 

Let k = conductivity, <r = specific heat per unit volume, 

A = cross-section of ring, c = its mean radius, 

s = an arc of ring measured from the heated spot, 
<j> = the angle subtended by s at the centre, so that s = c<jf>, 
H = heat supplied at the spot per second. 

The general equation of conduction is 

8 being the temperature. 

If the heat has been supplied for a long time t, the temperature will be 
sensibly uniform throughout and equal to 

m/27rcA<r. 

If we include the effect of conduction we may satisfy Equation (1) by adding 
to this a term of the form B<jb 3 ~f-C<£-f D, where 

B - H.c/27rAk 

and the other constants are undetermined. 

To complete the general solution we have to add terms of the form 

L + 2M £ cosi<£<r Ai (2), 
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We have also to satisfy the conditions that the flow of heat outwards at the 
origin is constantly equal to |H on either side, that is, corresponding to 
<£ = 0, cj> = 2tt. 

All these conditions are satisfied by the expression 

^=^^ + -^ 7 (0--7r) 3 + L + 2M,eos^6-^ (4). 

It will be noted that the last terms give no flow at 0, on either side, nor at 
the diametrically opposite point, as ought, from symmetry to be the case. 

The constants L, M are determined in the usual way by taking = 0, when 
t = 0, all over the ring. Thus 

r2rr TT/> /*2rr 

J o J o 

cos 2 icf> d$ + 7 -^ J j (<£ - irf cos i<j> cUf> = 0. 

J o 

And, finally, 

0- m + Hg ^"" 7r ) 3 ~- ^ Hc -^v 1 cos id> e-u (ft 

Case 2. Infinite Straight Bod Heated at the Ifiddle. — By putting 

c<p = x, i/c = 2?, 1/c = dz, and £ = oo , 
we pass, in Fourier's way, to the case of an infinite straight rod. 

n HX . HC7T 2 H r°°l -2.M/rr7 

= — - — -=+- — - — — r-r- ~ cos zx e~ ltz < a dx: 
2-rrAk SirAk irAh] z' z 

The second term, which is infinite, may be treated thus — 

tt 3 /6 = l/l 3 -f-l/2 3 + l/3 3 + . . . 



1 f°° 

c7r 3 /6 = ~ 2c 3 /^ 3 = I 2 3 <#£. 






If we integrate the last term by parts, dropping the factor H/ttAJc and 
include this term, we have 

f _ 1 + leos zx e-W'l °° + C XSmzx e- k ^dz + — f °° cos «c e~ u ^dz, 

L z z «L '■ s <r J 

o J o J o 

and the part outside the integral sign now vanishes at both limits. 
If we put n = Jc/cr, 



7T 



J n 2\A* 



1904.] Finding the Conductivity for Heat. 43 



and thus finally 



o z 2 ^ ni K 



#- -J + 4(2^fHJV^*)4=: ......... (6). 

2 Ah irAkl J J 2 s/ nt 

o 

Case 3. iJwiy with Radiating Surface. — If P be the perimeter of a cross 
section, and h the constant of radiation, the equation of the flow of heat along 
the ring is 

*A#0 _ A<r^?+PA0 (7). 

If the flow of heat from the source were continued for a long time t, the 
final temperature, neglecting conduction, would be given by 

2mhP v ' 

Including conduction, it is easy to see that the general Equation (7) is 
satisfied by an expression of the form 

z=z ~~J^6- /lP */ A *+B^ (8), 



where m 2 /c 2 = P/&/A&, 



1 (<)} 



The condition that the flow at either side of the origin should be equal to 
H/2 requires that 

— B . m/c . &A (<r w * — ^) = H/2, 

when t = 0, == everywhere. 

Thus 

L -— i_ + B (*» <*"*>+ «-* (+— >) + 2M ; cos ^ = 0. 

Determining the constants L r Mi in the usual way, 

H c 



L = 0, Mi 



7rA7H 2 -f-m 3 
Thus, finally, 

27rcAP 2 mAJc e mir —e- mn irAk i 2 + m 2 
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Case 4. Infinite Radiating Rod. — To pass to this case, put ccf> = x, c = co . 
it then w 2 = rh/Afe, m = c#, and — ■ = ■ — = e &* 

■*■ * 7 L1 wmir s>~mir 1 g—2cpjr 7 

f) — e ~~ PX — fi -mt/A<r[ COS ZX ntz 2i /-t-ix 

2Ak p ^Ak e - \ i /+f { * 

This integral may be expressed in terms of two integrals of the error- 
function type, but it is unnecessary to give the reduction here. 

Case 5. Thin SpJterical Shell.— The case of an infinite cylinder heated from 
a wire passing along its axis is the same as the two-dimensional problem of 
a plate heated at a point. 

To treat this, consider first the case of a thin spherical shell, heated steadily 
from a small circular hole, and take the boundary of the shell to be impermeable 
to heat. 

The general equation of the flow along the shell is 

h d \ /i n\ du ~] do /■% c\\ 

a^{ (1 - f " ) T^ = <r dt (12 >> 

a being the radius of the shell, /* = cos 3, $ = angular distance from hole. 
If e be its thickness, the final state, after a long time t, is given by 

Hi! 



*/ = 



4z7ra 2 et 



To satisfy Equation (12) generally, including this term, we require a term q 

such that 

H 



*i{«-<} 



CO 



dfil dfjbi 4:7re ' 

whence (I—/* 3 ) -=% = a (a-f C), where a = - — - . 

^ dfi ' 4zireh 

And the complete value of 6 will be 

CO 

o 

where X = i(i-\-l)kja?a, P$ being the usual zonal harmonic of degree i. 
The flow across any small circle $ is 

— he x 2ira sin # — — = 2irke (1 — u?) — . 

add- dfjb 

The part of this depending on P* vanishes when # = or fi =. 1 ; and, if 
C = 1, the part derived from q will be equal to H. Thus 

da a H i /1 \ 

4 = (w)' or *=-5a I °8< 1 -'* 
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The complete expression for is therefore 

iJL-iSfc log(i -^ + ? MiPi<rAt < 13 >" 

To determine Mj we have, when t == 0, 

+ 1 r+1 



M 



Pi%6-f] gVidfju = 0. 



Hence 

)^ + i(^i)fp^ = {(i_^)(p,|-j^)}^ 

The first term on the left hand vanishes except when i = 0. On the right 
hand (1— //, 2 ) -2 P { = a(l + /^) Pi = 2a, between the limits ; and 

(1-^) q **} = (!_/,) logCl-^x (1 + u) *&, 

which vanishes both when /*, = 1 and when /i = — 1. Thus 

Mj=-o(2»+l)/i(t + l> 
The constant term M is given from 

M Fd/i = a log(l-/*)d/* = — «{(l-/»)log(l-/»)-(l-/*)}iJ 

= a (2 log 2— 2), 
or M = a (log 2 — 1). 

Thus the complete expression for is 

where a = H/^ireh 

When t = 0, 8 = everywhere, and this requires that 

lo g lz^ + i + s^±ip i = 0> 

B 2 i(i + l) 

an expression which may be readily verified from the expression 

(1 - 2 fix + a?»)-* = 1 + 2P^. 



46 Prof. C. Niven. On a Method of [Dec. -5, 

Case 6. — Infinite Plane Plate. This case, which also includes that of an 
infinite cylinder heated along its axis, may be deduced from the foregoing 

one by putting 

' r \ 3 



1— /£,= 2 1— I 



r being the distance of any point from the centre, and by taking 

ija = z, 1/a = dz, \ = kz 2 /a ; 



also Pi = J 


o O0> 




and the expression for 9 is 




0= H i 7 (21o g ; + l) + 9 Hl 7 - 
4firk \ 2 / 2?r& 


log a — 


r 00 "l 
st 1 ^ (rs) e~ ntz \lz 



(15), 







where B^ stands for the heat supplied per unit length in unit time, and 
n = k/cr. 

The integral, which itself is infinite, is coupled with an infinite term, the 
difference between the two being the finite part on which the solution 
depends. The same point came up in the deduction of the solution for the 
infinite rod from that of the ring, but in that case it was possible to assign a 
finite form for the difference. 

The form of the integral in (15) suggests, as will be seen later, another way 
of obtaining the solution, but it is possible to deduce from it a form of 
solution satisfying all the conditions of the problem. For if we differentiate 
the above expression with regard to the time, we have 

rW H r°° 

or, putting f = zs/nt, x = rjy/nt, 

we have -=- = -r-V T u > w = ?J (#£) e-?d£. 

at 4:7r/ct J 

It is easy to verify by direct differentiation that 

Z—+UX = (J, 

ax 
and therefore u = Ae~ x ^ 4 

To find the constant A, we observe that when 



x = 0, u = I £e-&d% = | ; so that A = \ 



o 
Thus ^ = ^i 1 _ :c2/4 e = H^ ? r \ e -^ dt 

dt kirkt 4nrk) 



o 
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But a? = r^jnt, and therefore 2dx/%-{-dt/t = 0. 



Thus 



tij 



o 



-00 



Another Solution of this Case. — The expression just given for 6 was 
originally obtained by a different method, which may be now given. The 
integral contained in the Equation (15), viz., 



sr 1 J (rz) e- ntz2 dz f 

J A 





may be written, using the same symbol £ as above, in the form 



?- 1 J (fr/ v^5) e-Pd£ 





which shows that 6 is of the form <f> (r 2 /nt). 

But the general equation of the flow of heat in two dimensions is 

m d?e = ide 

dot? dy% ndt 
which may be satisfied by = (f>(v), where v = (x^ + y^/nt, if 

A dd> t 4 d?<b dd> 

4- r + 4?;-— x = — v-P- , 

cm? a^ at? 

and therefore ?; ~- = Ae~ v ^ t 

aw 

The total flow across the boundary of the circle of radius, r, in 

= — 2irrk -e- = — 4tt^ -=£ , 
&r aw 

and when r = 0, or v = 0, this is to be equal to H 1 . Therefore 

A == —Hij f 4:7rfo 

TT -.00 

Thus 6 = -pi. I r^-^, 

47r/£ 
*'o 

the upper limit of integration being taken to make 6 everywhere = 0, when 
t = 

Writing # 3 for v, this expression becomes the same as that previously found 
in Equation (16). 
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Case 7. Case of an Infinite Solid with small Spherical Source. — The equation 
of motion in this case is 

\clr 2 r dr) dt 

d*(0r) d(6r) „,-, 

where n = 7^/cr. 

When t is infinite and the flow has become steady, 

= A/r + B, where h . 4ttA = H. 
But the general Equation (17) is satisfied also by 

0r = (z) cos £r e~ ntz *dz } and by <£ (2;) sin 0?' e" ntz2 dz. 
**o Jo 

If we choose the latter, and put <£ (z) = C/», we may take as the solution 

applicable to the present case 

J o 
For, 



(1) When £ = 00 , it reproduces the steady state ; 

(2) When t = 0,0 =±(^ +C f !i nOT 



cfe 



2 




r(« + 2 /' 
which vanishes if C = — 2/tt .H/^ttJc; 

«r 7T I 7T J 

J 

which is equal to H when r = 0. All the conditions are therefore satisfied. 
The solution therefore is 

JLJL Jl I -* £i J feixl &'l m-tni 7 \ XJL JL / -I JL 



= 



o v 



